“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1972-06 


The wavelength dependence of the 
modulation transfer function of the atmosphere. 


Mackin, Jere Gene 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/16436 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


THE WAVELENGTH DEPENDENCE OF THE MODULATION 
TRANSFER FUNCTION OF THE ATMOSPHERE 


Jere Gene Mackin 








NAVAL POSTGRADUATE StHOOL 


Monterey, California 





THE WAVELENGTH DEPENDENCE 
OF THE 
MODULATION TRANSFER FUNCTION 
OPStHE ARMOSPHERE 


by 


Jere Gene Mackin 






Thesis Advisor: S. H. Kalmbach | 


per ey T148516 


Approved for public nelease; distribution unlouncted. 





The Wavelength Dependence 
of the 
Modulation Transfer Function 


of the Atmosphere 


by 


Jere Gene Mackin 
Lieutenant Commander, United States Navy 
B.S., United States Naval Academy, 1964 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN PHYSICS 


from. che 


NAVAL POSTGRADUATE SCHOOL 
wine “1972 





ABSTRACT 


An apparatus was designed and constructed to measure the modulation 
transfer function of the atmospheric medium as a function of wavelength 
of light. The modulation transfer function (MTF) in this case was 
considered to be the ratio of the modulation of an optical target at 
a fixed range to the modulation at a range short enough so that atmo- 
Spheric attenuation was negligible. 

Modulation measurements were made using a telescope which focused 
the image of the target on a slit. The image was mechanically scanned 
across the slit by a rotating accor and the intensity at the slit was 
measured by a photodetector. 

The preliminary results indicated a linear decrease in modulation 


oO 
of 0.1 per 1OOOA. 
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I. INTRODUCTION 


With the increasing complexity of optical systems in recent years,. 
it has become necessary to find a general figure of merit by which to 
compare them. By viewing the object of sacs as an cone and the 
image formed as an hue, it has become common practice to use the 
theory of the transfer function as the descriptive link between the 
two. 

Generally speaking, the quality of an image formed by an optical 
system is a function of three things: 1. The total abberation in 
the ray optics design of the system, 2. The (inherent) diffraction 
caused by the wave nature of light, and 3. Manufacturing inaccuracies. 
Paeetewedse where inanufacturing and abberation effects can be consider-— 
ed negligible, the system is said to be diffraction limited. The dif- 
fraction limited case is the cloSest to ideal that can be reached. 

The observable effect of diffraction is a reduction in contrast which 
increases with the spatial frequency of the target. 

When dealing with the intensity distribution of a periodic target, 
an important comparative parameter is the modulation, defined as the 
difference between the maximum and minimum —— amplitude divided by 
the sum of the two. This parameter is the ratio of (twice) the ampli- 
tude variation to (twice) the average intensity. 

In the most common usage, the plot of modulation versus the spatial 
frequency of the target (cycles per distance on the target) is anes 
the modulation transfer function (MTF) curve. Since this is neared 
to 1.0 at zero cycles, it will be defined in this paper as the ratio 


+4 


of the modulation at a fixed distance to the modulation at approximately 





zero distance for a fixed spatial frequency at a particular wavelength 
of received light. Since the atmosphere is the common medium of 
propagation for most optical systems, it is necessary for both pre—- 
diction and evaluation of system performance that an understanding of 


its characteristics be applied to MTF calculations. 


Poemeeoky OF Hit FOR A DIFFRACTION LIMITED SYSTEM 

The response of an optical system is determined for basic patterns,, 
such as points, eee or edges. The responses to these patterns are 
called the point, line and edge spread functions respectively. For 
responses to complex objects, such as periodic targets, the optical 
transfer function is used if both amplitude and phase are of concern, 
and the MTF if amplitude only is desired. Since any object may be 
Cem@erdered a collection of points, there is a connection between the 
point spread function cual the transfer functions. 

The point spread function is’ the normalized diffraction pattern 
formed by a point object. It is usually plotted as illuminance as a 
function of distance in the image plane. 

Once the point spread functions for a sensor system are known, the 
illuminance of any object can be obtained by summing the point spread 
functions of the target object. This process can be represented by 
a convolution integral. As a general example, consider an object plane 
(x, ) and an image plane (x, ). Let the function O(x, ) and the spread 
function SCx,) define an object. ‘The illuminance of the image, I(x); 


is given by the convolution of the object functions: 


oo 


I(x,) = § 862,)0(x-%,)dx, 


Lo @] 








Figure 1. IJilustration of convolution process. 


Figure 1 shows the convolution of an object function and a spread 
function to produce an ee function. The image function does not 
exactly reproduce the object because only three points in the image 
Function are illustrated. The appearance of the argument (x, -x,) 
implies that one function is translated or scanned along the other 
Function so that each point in the object is traced through the 
system to the image, and the spread function is then scanned over 
the image point. 

For a linear system, the response due to a number of inputs is 
the product of the transfer functions of each. In an optical system, 


thie mastematical feundarwen Ws the Pourler™transforn: 


G(w) = { elsdexp(-iwxyax 
In any real system, the spatial varlable is usually at least two- 
dimensional. Taking the Fourier transform of both sides of the con- 


volution integral, 





I(w) if S(x, )0(x, -x, )exp(-iwx, )dx,dx, 


f Sx, )dx, f O(x, —x, )exp(-iwx, )dx, 


-—O 


J stx)ex, foCarexp(-iadexp(-iwe,) du 


—-& ~—~c 


where u = x_-X 


= c 2 
= J 8(x,dexp(-in%, 2x, focwexp(—iwa) du 
So that I(w) = S(w)0(w). 
The spread function can then be written $(w) =a When. normalized,,. 


the Fourier transform of the spread function is the optical transfer 


co 


PMC TION . 


S(x)exp(-iwx) dx 
OTF(w) a 


~S(x)dx 


=-CO 


The MTF as a function of k = oe can be obtained by: 


MTF(k) = lore) : 


It should be pointed out here that k is the spatial frequency in lines 
an 

per distance, and that w is therefore a spatial angular frequency, in 

radians per distance. Recalling that 


exp(-iwx) = exp(-i2nkx) = cos(2nkx)-isin(27kx), 


oO oO 


| ( s(x)cos(2nkx)ax-i f S(x)sin(27kx ) dx 


MnO) <i> | ae ; 
f S(x) dx | 


The numerator is the vectorial sum of the cosine and sine Fourier 


transforms of S(x). 





As an example of the calculation of the MIF of a spread function, 


Gensider a function defined as: 


-a a 
ae a) 
S(x) = 
0, elsewhere 


This might represent the spread functon of a scanning spot of diameter 


a with spatial frequency =. Since the spread function is symmetrical,. 
a/? 
f Saint 2meodx = 0 for all.k 
~a/2 2 
/2 | ‘~ 
S(x)cos(27Tkx)dx = 2 f cos( 27kx) dx 
-a/2 
~a/2 / ay? 
- 9 sin( 27kx) 
7 Onk 
0 
_ sin(tka) 
mk 
a7 2 


The denominator, f S@ dx = a, so that 
-a/2 


MTF(k) = | sin(tka) 


tka 





It can be seen that the MTF vanishes at integer multiples of =. 

The existence of the MIF beyond its cutoff frequency of = is 
referred to as spurious or false resolution. This can be explained 
by considering a bar chart of lines and spaces of width a. If the 
scanning spot of diameter a is positioned so that the center of the 
spot traces the edge of each successive line, the output of each 
scan will be the same because a will scan the same amount of line 
and space. If the spatial frequency is increased, the lines and 
Spaces become narrower, and therefore closer together. In this 


case, the amount of line and space covered each scan varies from 





nearly all line to the opposite, and the varying output looks like 
resolution, although no useful information can be obtained from it.. 
In the case of a diffraction limited lens, the spread function 
is the pattern consisting of the central Airy disk surrounded by 
pRegwessively fainter rings. The MIF is the Fourder transform of 
this diffraction pattern. Consider the general case of Fraunhofer 
diffraction of a circular aperture. Figure 2 illustrates the co- 


ordinates used. 


Figure 2. Fraunhofer diffraction of a circular aperture. 


It is desired to find the illumination at a point P! due to a 
source at point P after the radiation has passed through an aperture 
A. The amplitude of P' is found by summing all the contributions 


from the area elements in A. 


a . t e+ i 
Ute!) = otf uCKY) epi dA = CC! f explites ) )dA 


Because kr is large, the numerator changes rapidly with r, while the 
denominator changes slowly over the range of integration, so long 


as rr' >> Pi ee , 


18 





ly 2 
[(x-x)2 + (y-¥)? + 22] 


r= 
, 5 1/2 
i -2xX - 2yY 
O 2 
z 
O 


ee ay et = Ot y¥)* 
2 P Fang or 3 
O 
A similar expression can be obtained for r‘'. In these expressions,. nr. 
and rs are the distances from the origin to P and P' respectively.. Now 


let 1 and m OKs TENE Gest tivo @lrection cosines for Peand L' andam' be 


those for P'. These cosines may be written: 
ieee ae) , ae 
eee eee 5 tT Re yr 
O O O O 
Substituting, 


U(P') = Cexp(ik(r, + ep ff exp fikC(1 —- 1')X + (m-m')Y 
; Ae 


Z 2 
ty t 
oy Ly cee sy?) - eT axa 
2 a dl 


r ee 
O ‘e} 


2 Z 
= see ax, the terms in the second 


t 
© 


In the case where r 


3 














r 
O 
bracket may be neglected. This is the Fraunhofer diffraction approxi- 


mation. The amplitude expression may be simplified to: 


U(Al,Am) = off expLik(XAl + YAm) JdxdyY 
A 
Assuming the aperture does not transmit perfectly, and calling its 


transmittance T(X,Y), 


Dene c ff T(X,Y)exp(ik(XAl + YAm))dXd¥ 
The amplitude of the diffracted wave is the two-dimensional Fourier 
transform of the aperture transmittance, provided the aperture co- 
d 


. ' al ; 
ordinates are measured in units of ee Using the Fraunhofer 


aia 





approximation, the linear coordinates may be written, x = rial 


_ t 
y = ram 
Tee 7 TX, YexpGk (ax + y¥))axay 
roe) O 


Consider a circular aperture of radius a. In arder to integrate 


in polar coordinates, we substitute 





X = Reosé, Y = Rsiné@ 
xX = rcos®d, y = rsing 
dA = RdRd@, ff = rt 

J ae 

ikrR ; z 
e=e J R f exp | = (cos @cosd + sindsing) | dédR 
O O 
a 27 





cos(@-¢) dé@dR 


S 
i 
©) 
2 
on 
oO 
us 
Ae) 
aN 
=e 
Khl a 
ne) 


27 
Recall that f exp(ixcoSg)da = 2nd (x), where J 6x) is the zero order 
O : 


Bessel function. Since the integrand is periodic, we may add any 
finite amount to both limits simultaneously without changing the 

- krR 
value, so that U = ave f RI (AS )ar. The recurrence relation for 


pesese) tfumetions states that 


d gaps eae tl. 
ae Je? =a JO) » SO that 
fa (ides = x (x) 
O ai 
oJ (Ka) 
5 ee oe 
Gos (qe se), == 
kra 
a 


The irradiance is the square of the amplitude, or 


i 








mC 
aes AP 
etre kra I_, where I = 


Ph =sdpemeure area, — = Local power 


The corresponding distribution 1s called an Airy pattern, and the 
area enclosed by the first minimum is an Airy disk. The diameter of 


the disk turns out to be approximately 1.22 ee » where a’ is the 


na 


apex half angle of the cone converging on the image point. 


To determine the amplitude spread function, consider Figure 3.. 





projecties of Y ow 


peference yphare 


Figure 3. Derivation of lens spread function. 


ae 





The spherical wavefront of radius R has its vertex, A, at the 
center of an aperture, Q, which can be considered the exit pupil of 
a diffraction limited lens. pcan) is a point in the image plane.. 
The amplitude of the wave at P is obtained by integrating from all the 
area elements dx'dy' over the wavefront. For-each point B'(x',y') in 
the ee a projection is made to a point B(x',y') on the wave— 
front. A wavefront area ds corresponds to _— dx’dy'. The projeetion 


1s made such that cos?6. The amplitude Uy at P can be 


eee. = 
qe ay’ 
obtained by integrating over the area elements at Q. 


~' 


1/2 
Uy) =f [Eexp(iks)cos36 dx'dy'] 
$ M2 7" 
= a oxpc nani expl (=) (xx' + yy")] dx"dy' 
E = aperture illumination 


@ = total energy flux in the aperture 


50 = aperture area 
: alg : : E : 
Assumptions: RR pees = entre 1m nation. 


The integral above is the Fourier transform of the aperture function. 


2g, (krA) On 
For a circular aperture this becomes U (r) = ————— , where k = —, 
O krA ro 
4 is the wavelength in vacuum. A is.a numerical representation of the 
172 
aperture, and r = (x* + y%) . Note that this is the square root of 


the Airy function. A pupil function can be defined, V(x',y') = 
T(x' ,y' Jexp(ikv(x',y')), where T is the amplitude and v the phase 
deviation in the aperture. Writing the amplitude spread function in 


terms of the pupil function, 


co 


f V(x" y" JexplSS x" + yy')Jdx'dy' 
2 i < -_ - 


f VG jas" dy. 


14 


U(x,y) = 





The amplitude spread function is, then, the inverse Fourier 
transform of the pupil function. 

To obtain the transfer function of a lens, it is necessary to 
take the Fourier transform of the spread function. New variables 


t t 


i Lt =, ae = 3 
are defined; ee ip? and Ss ip > SO that 


= t ' t & ' t 
UG..,) = c ff V eo re iCxv} + yvy dv dv, 
-42p2 


ee 
ff Neo align dy 


This is the inverse Fourier transform of V', the pupil function,. 


Cy. 


written in the new variables. For an incoherent object, the transfer 


function is the Fourier transform of the intensity spread function,,. 


OLE(x.y) . rE * 
M0%5) SETEGG YT mG yy) . where E(x,y) = lucy) 


< a vf —— _] ds 
B@ey) = Ciee (Vv) ih (v")] 
Applying the convolution theorem of the Fourier transform, 


CR pe Cel Hh hyd} 
x? “y # led 


ed Oe pare = ry ye 
Sc Ra 


[EI 


See (VG@)1}" = Vv (-x) 


Thus the convolution changes to correlation, 


co 


Gee See ene £(x)Oa(X) =f rice ent 


x? -y di (FE) 


The denominator, di (E) = Viov'= (0,0) = "OT, “fiewtotal Pitk through 


: Vovi= , 
the aperture. When normalized, pe ea = On, which reépyvescemec 
i 


aS 





the value of oe integrated over the common area when two figures 


of identical cross section are displaced by ens Taking the example 


Conmion aveca 


i i lar apertures a h Je SS, 
Of Preere 4, where two eircular apertures are shown, the MIF = 





Figure 4. Correlation of two identical circular apertures. 


v v Vv 
From Figure 4 it can be seen that rcos? = eae so that ceso = => a 
2 
ae - (rsind)(reos¢)] = 2r*($-sindcosd) 





The common area is given by 2[ 


MTF = 


2 ($-singcoss) 


If the maximum spatial frequency is Ko? then the normalized spatial 


k Vy 5 1/2 
frequency = — = — = cos, and sind is equal to (l-cos“¢) . 


“lk k eae 2 
[ cos Cae 6 Cc? j 
O O O 


MTP = 


23, (tk. x) 


This is the Fourier transform of , the Airy pattern, so 


WK OX 
that the MTF is that of a diffraction limited lens. The interesting 


point is that the MTF is the convolution of the aperture with itself. 


B. THEORY OF THE EXPERIMENT 

The wavelength dependence of the MTF of the atmosphere was to be 
obtained by optically scanning a 4' x 6' bar chart made of alternating 
black and white strips. Each strip was ten centimeters wide, with a 


resulting spatial frequency of five lines per meter. The angular scun 


Als 





provided by a rotating mirror was a total of 0.0038 radians, and was 
interrupted by a mechanical chopper at the rate of 29 chops per second 
to provide a reference Zero. The chopper and rotating mirror were 
constructed as shown in Figure 5. The wavelength of the received light 
was to be varied through the visible and into the insiemericiae and the 
modulation observed and referenced to the modulation obtained at 


essentially zero distance to yield the MTF. 


C. SVOLUTLON OF THE EXPERIMENTAL APPARATUS 

The block dneteeen of Figure 6 was a common factor throughout the 
Woseemen effort, Meer several modifications were made to the optics in 
an effort to improve signal amplitude and accessibility of the slit for 
visual target monitoring. 

iitemicthira! Optical arrangement was as show in Figure 7, using 
the monochromator Seren of a Perkin-Elmer Model 13 infrared spectro- 
photometer. In this arrangement, the light passed through the entrance 
slit, was dispersed by the prism, passed through an array of mirrors 
and slits, and was finally reflected up through an aperture in the 
cover of the monochromator to a silicon photocell by mirror M6. The 
detector circuit was as shown schematically in Figure 8, except that 
the rectifier diode had not yet been added. Wavelength scanning was 
achieved by changing the angle of M2 to cover the entire prism 
dispersion pattern. 

A problem was discovered with this arrangement in that the losses 
suffered by the light in passing through 9 mirrors (including two in 
the telescope), through the prism twice, and into the detector were 
SO great that the signal amplitude out of the detector was too low to 


be observed above the noise on the Tektronix 451 oscilloscope. The 
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first effort to improve this situation was to shift from a Type E 

Low Level Dyisciearconyaieel AC Preamp plug-in scope module (which had a 
lower range of 50 microvolts per cm) to a type 1A7 High Gain Dif-— 

ferential Amplifier (with a lower range of 10 microvolts per cm).: 

This was still unsuccessful. With the entrance slit full open at 

2 _— the white bars on the target were visible on the scope,. but 

the black bars were lost in the noise. 

The next modification was to bypass the monochromator entirely,,. 
and place the detector directly behind the entrance slit.. This gave 
detectable signal amplitude, but a problem of overshaot on the scope 
coupled with capacitive discharge through the detector circuit. 
produced a negative response of approximately the same amplitude as 
the signal. At this point, it was decided that at least part of the 
problem could be solved by dealing with the capacitor. The value of 
the capacitor was changed twice, ending with .047 microfarads.. It was 
also observed at this time that the photocell output as observed on 
the scope was overshooting in the negative sense on. each chop, either 
eeea result of the internal escilleseepe er detector circuit charac— 
teristics. 

By partial dismantling of the system, it was possible to determine 
visually that the chop on each cycle was occurring while the image 
of the target had been swept past the slit, which implied that there 
was only background present when the chop occurred. This problem was 
corrected by rotating the cam on its shaft. 

At this point in the research, an attempt was made to change 
detectors from the silicon photocell to an InSb photovoltaic, liquid 


nitrogen cooled detector manufactured by Barnes Engineering. In orde? 


22 





to accommodate the detector, shown in Figure 9, the optics were 
rearranged as shown in Figure 10. The documents pertaining to 

the detector failed to show the associated circuitry required, and,. 
as a ne. no useful data were obtained from it. An attempt was 
made to take the output directly from the detector to the scope, 
Since it operated at zero volt bias and had low maximum current and 
voltage ratings, but no signal was aurea: 

In order to alleviate the capacitor discharge problem, the 
rectifier and resistor were added to the detector circuit so that 
any discharge would occur through the internal circuit vice through 
the output. 

Since phone calls and local efforts were unsuccessful in getting 
results from the InSb detector, the silicon photocell was reinstalled 
using the optics as shown in Figure 10. Wavelength variation was 
achieved through the use of interference filters. The results obtained 


will be presented in Section II of this paper. 
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II. RESULTS OF THE EXPERIMENT 


The following data were obtained using a target illumination of 
approximately 200 watts per square meter by incandescent jigs a 
chopping frequency of 29 per eoeasar and essentially zero attenuation 
due to distance. 

The reference modulation of the target based on 46 measurements 
at distances between 20 and 120 meters was 0.647 + .011. This should 
be normalized to 1.0 in dealing with future data. Graph 1 is a plot. 
of modulation versus wavelength at zero range. The reliability of 
this graph is questionable since the interference filters used were 
of varying transmittance. In order to resolve this problem, a set 
of matched, high transmittance filters should be used. It ls expect— 
ed that the trend of this graph will reverse to some extent at the 
higher wavelengths, using an infrared detector at greater range, due 
to the infrared window of the atmosphere. 


Figure 11 shows scope: photographs of typical signal presentations. 
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Figure lla. Typical 
oscilloscope presentation. 
Highest amplitude is white 
Bar, double peaks are black, 
and chop is between. 
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Figure llb. Expanded 
BEGOpe Trace. 
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